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ABSTRACT 

We present an analysis of recent XMM-Newton and HST archive data of the 
ultraluminous X-ray source NGC 1313 X-2. Quasi-simultaneous observations 
taken with XMM-Newton, HST and VLT allow us to study both the X-ray light 
curve and its correlation with the optical emission of the two proposed ULX 
counterparts. At the end of December 2003 the source experienced a short, but 
intense flare, reaching a maximum luminosity of ~ 10^° erg/s. At the same time, 
the optical flux of both the suggested counterparts did not show pronounced vari- 
ations (< 30%). Assuming that the ULX emission is isotropic and taking X-ray 
reprocessing into account, the optical data for one of the proposed counterparts 
are consistent with it being an early type, main sequence star of ~ 10 — ISM© 
losing matter through Roche-lobe overflow onto a ~ 120Mq black hole at an 
orbital separation corresponding to a period of ~ 2 days. 

Subject headings: galaxies: individual (NGC 1313) — stars: individual (NGC 
1313 X-2) — X-rays: binaries — X-rays: galaxies 



1. Introduction 

When, at the beginning of the 1980s, po int-like, off-nuc lear X-ray sources in the field of 



nearby galaxies were first detected (see, e.g.. lFabbiandll989l ). it was immediately recognized 



that their luminosity was unusually large. If physically associated with their host galaxies. 
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these sources would have an isotropic luminosity in excess of the Eddington limit for a IOMq 



Roberts & Warwick 


2000; 


Colbert & Ptak 


2002; 


Swartz et al. 


2004; 



l2005h 



It is estimated that a s ignificant fract i on of ULXs are interacting supernovae or back 



ground AGNs (~ 50%; see iFoschini et all l2002bl : iMasetti et all l2003l : ISwartz et al.ll2004[ ). 
However, the X-ray vari ability of many of th e m is similar to that observed in Galatic 



X-ray binaries (see e.g. iLa Parola et al.l 12001 



2001 




Colbert & Ptak 


2002 




Swartz et al. 


2004 



Zampieri et al.ll2004l . hereafter Z04). The recent detection of a 62 days modulation in the 
light curve of M 82 X-1, interpreted as the orbital period of the system, provided a dir ect 
confirmation of the binary nature of at least some ULXs JKaaret. Simet fc Lang|[2006al lb[). 
Moreover, ULX spectral properties share simila rities with those of Galactic black hole X-ray 



binaries (BHXRBs, e.g. iFoschini et al.l l2002al ). In several cases the spectrum can be well 



reproduced by a multicolor disk (MCD) blackbody plus a power law (PL), although the tem- 
perature of the MCD component is often mu ch lower than that observed in BHXRBs (e.g. 
Miller et al.ll2003l . l2004j : IPeng fc Kaaretlboosi ). For the brightest ULXs, a possible curvature 
above 2-3 keV has been recently reported and more sophistica ted spectral models appear to 
give better agreement with observations (IStobbart et al.ll2006l ). 



All t hese properties along, in some cases, with the detect i on of stellar-like optical coun- 



terparts ([Roberts et al.ll200ll:lGoad e t ahlbooi Eiu et allbood I2OO4I: 



Kaaret 2005; Mucciarelli et al. 



2005 



Kaaret et al. 2004; Z04 



hereafter M05; ISoria et al. I l2005i r strongly suggest that 



a sizeable fraction of ULXs are accreting X-ray binaries. The present debate is focussed on 
understanding what type of binar ies they are . Many of the ULX properties can be ex plained 
if they do not emit isotropically ( King et al. 2001 ; King 2002 ; King fc Pounds 12003) or are 



dominated by emission from a relativistic jet (IKording et al.l l2002l : iGeorganopoulos et al. 



2OO2I ; iKaaret et al.l 120031 ). Another possibility is that they are truly emitting above the Ed 



dington limit for lOM^ , either because accretion proceeds through a slim disk (lEbisawa et al. 



2OO3I ; iKawaguchil |2003| ) or because the comp act object is an Interniediat e Mass Black Hole 
(IMBH) with a mas s in ex cess of IOOMq (e.g 



Colbert & Mushotzkv 


1999; 


Miller et al. 


2003 



Patruno et al.ll2005l . l2006l ). Despite the inherent difficulties due to the low counting statis- 
tics. X-ray timing analysis has been attempted in some ULXs and led to the detection of a 
quasi periodic oscillation in the power density spectrum of M 8 2 X-1. This may represent 
a powerful, independent rnethod to measure the black hole mass (jStrohniaver fc Mushotzky 
2OO3I : IPiorito fc Titarchuklliooi : lOewangan et al.lbood : iMucciareUi et allbooeh . 



Multiwavelength observations are an unvaluable tool to investigate the nature of ULXs. 
Radio emission, when present, gives important clues a bout the geometry, energetics and 
lifetime of ULXs ( iKaaret et al.l |2003| ; iMiller et al.l |2005|) . Optical observations are crucial 



3 



to identify ULX counterparts and to study the properties of putative ULX binary sys- 
tems. Up to now only a very small number of ULXs have been coiivincinKly associated 
with stellar objects of known spectral type (e.g. iLiu et all 120021 . 12004 : iKaaret et al.ll2004j : 
M05). All these ULXs are hosted in young stellar environments or star-forming regions 
and their optical counterparts have properties consistent with those of massive stars. Some 



ULXs are also associated w ith extended optical emission nebulae (jPakuU fc Mirionil 12002 
Pakull. Grise fc Motch|[2006[ ). 



NGC 1313 X-2 is a well studied ULX. The X-ray variability, high (isotropic) luminosity 
and presence of a soft X-ray spectr al component make it a prototypical object. Furthermore, 
the presence of an emission nebula (jPakuU fc Mirionill2002l : iRamsey et al.ll2006l : Z04) and the 
detection of optical counterpart (s) (Z04; M05) provide a considerable amount of information 
on the ULX environment, available only for a very limited number of objects. Here we present 
a systematic study of the X-ray and optical variability of NGC 1313 X-2 based on archive 
data of the XMM-Newton satellite, and the ESO VLT and i/S'T telescopes. Observations are 
reported in § |2] and results in § |3l where also a model for the optical emission is presented. 
Discussion follows. 



2. Observations 

2.1. X-ray observations 

XMM-Newton observed NGC 1313 in 9 exposures taken between 2000 and 2004. All the 
observations are listed in Table [H A detail ed analysis of the O ctober 2000 data of NGC 1313 
X-2 have been performed by Z04 (see also TuroUa et al.|[2006l ). Here we report results from 



an analysis of the EPIC pn exposures of the 8 more recent observations. Data reduction 
and extraction have been carried out with standard software (XMM-SAS v. 6.0.0). All the 
observations are affected by solar flares. The good time intervals left after subtraction of the 
high background periods (when the total off-source count rate above 10 keV is > 1.0 counts 
s~^ for EPIC pn) are reported in Table [H For the analysis we consider all the exposures 
with a good time interval longer than 1 ks. After performing standard cleaning of the event 
lists, we extracted source counts from a circle of 40" centered on the position of NGC 1313 
X-2 (Z04). The background counts were extracted from a circle of 50" on the same CCD. 

The spectral analysis was carried out within XSPEC (v. 11.2.0). A two-component 
model consisting of an absorbed multicolor disk blackbody (MCD) plus a power law (PL) has 
been employed throughout, as it is routinely done for BHXRBs. Pre yious applications of the 



same spectral model to a number of ULXs, including NGC 1313 X-2 flMiller et al.ll2003l . 12004 
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Cropper et al.ll2004j : iKong et al.ll2005l : Z04) gave a satisfactory fit to the data. The absorbing 



column density inferred from the different datasets is consistent with a constant value. We 
then performed again the fits fixing to the average value weighted by the exposure time 
{Nh = 4.02 X lO^^cm"^). For consistency, we also repeated the analysis of the 2000 EPIC pn 
spectrum of NGC 1313 X-2 following the procedure described above. The results from the 
spectral fits are reported in Table [2l The XMM-Newton fiuxes were consistently derived from 
the parameters of the spectral fits and are reported, with the corre sponding lu minosities, in 
Table [2] (a distance of 3.7 Mpc was assumed for the host galaxy; lTullylll988l ). The errors 
have been estimated from the maximum and minimum values of the fiux, obtained varying 
the fit parameters systematically. 



2.2. Optical observations 

B, V and R images and spectra of NGC 1313 X-2 were taken with VLT+FORSl in 
December 2003. The results were presented by M05. HST images of this field were also 
obtained with ACS in two epochs (see Figured]). The optical observations were performed 
in parallel with the XMM-Newton pointings. The observation log of the VLT and HST 
images considered here are reported in Table [T] Aperture photometry was performed on 
the drizzled cali brated data (reduced by the HST pipeline) and transformed to the Cousins 



system following ISirianni et al.l (120051 ) (see Table [3]) 



The uncertainty on the if5T magnitudes are dominated by the calibration error (~ 0.03 
mag for the i75'T photometry of point sources), including filter transformation. As a further 
check of the internal consistency of the HST photometry, we compared the magnitudes of 
thirteen field stars obtained in the two epochs. The difference is significant only for one 
source in the sample (~ 0.3 mag). Excluding this source, the variability of which is probably 
intrinsic, the magnitude changes are randomly scattered around zero, with a mean absolute 
deviation of 0.04 mag. 

The HST images clearly confirm that two distinct objects are present inside the X-ray 
error box of NGC 1313 X-2 (see Figure [1]), as first shown by M05 on the basi s of VLT data. 



If we assume Ay — 0.3 (Z04), taking Galactic absorption into account (the ICardelli et al. 



19891 extinction law with Ry = Ay/E{B — V) =3.1 has been adopted throughout), the 
unreddened colors inferred from the 1st HST epoch are (B-V)o ~ —0.13 and ~1.46 for CI 
and C2, respectively (see Table [3]). The color for object CI is consistent with that derived 
from VLT data ((B-V)o ~ —0.2; M05). On the basis of the HST photometry, the color 
of object C2 is close to that of a K3-K4 supergiant. For both objects, there is evidence 
of variability in the V band between the two HST epochs (~ 0.1 mag. Figure [21 see also 
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Ramsey et al.ll2006l ) 



3. X-ray and optical variability 

3.1. X-ray light curve 

Figure [3] shows the (unabsorbed) X-ray flux for all the available observations of NGC 
1313 X-2. The XMM-Newton data were derived from the best fltting spectral models (see 
Table [2]), while the Einstein, RO SAT and ASCA data are taken from Z04. 

Until 2000 NGC 1313 X-2 exhibited variability up to a factor of two on a timescale of 
months, with a maximum luminosity of ~ 4 x 10^^ erg s~^ (Z04). Around December 25, 
2003 (observation 6 in Tabled]), the source produced an intense flare, reaching a maximum 
unabsorbed flux of ~ 10.0 x lO'^^erg cm'^ s"^ (Table [2]). At the distance of NGC 1313 this 
corresponds to an intrinsic luminosity of ~ lO'^^erg s~^. Clearly this value depends on the 
adopted spectral model and hence should be taken with car e (see e.g. the slight y smaller 



values recently obtained for the October 2000 observation by lStobbart et al.ll2006l adopting 
a more sophisticated spectral model). We also measured the fluxes of another ULX in the 
fleld (NGC 1313 X-3), known to be an interacting supernova (SN 1978K), in order to check 
if the signiflcant luminosity increase was real or artiflcially produced by residual systematic 
effects between the two observations. The flux of the supernova (~ 8.2 x 10~^^erg cm~^ s~^) 
is consistent with a constant, within the uncertainties (the variation is < 20%). Hence we 
conclude that the luminosity increase of NGC 1313 X-2 is signiflcant and fully qualifles NGC 
1313 X-2 as a bright ULX. 

From the observed maximum luminosity [Lmax ~ 1.5 x 10^° erg s~^) and assuming 
isotropic emission, the black hole mass Mbh obtained setting L^ax = -^Bdd {LEdd is Edding- 
ton luminosity) is ~ I2OM0, about 2 times larger than that previously estimated by Z04. 
Sub-Eddington accretion would imply an even larger mass. 



3.2. X-ray spectral changes 

The spectra of numerous ULXs, including NGC 1313 X-2, are well described by a 
MCD-I-PL model, similarly to those observed in Galactic BHXRBs. In fact, there is also 
some evidence of a closer similarity, inasmu ch as some ULXs appear to show state transitions 



(IMakishima et al.ll2004j : IWinter et al.ll2005l ). In the following we summarize the main results 



of an analysis of the X-ray spectral variability of NGC 1313 X-2. The most signiflcant 
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result is that the slope of the PL component seems to correlate with the flux, i.e. at higher 
fluxes the spectrum hardens (see Table [2]). This behavior was already noticed by Z04 on 
the basis of a comparison between two ASCA observations and is opposite to that usually 
shown by Galactic BHXRBs. A sirn il ar co rrelation was also observed in a few ULXs in the 
Antennae galaxy by iFabbiano et al.l (120031 ). The MCD component is important in the Oct 
2000 and in the 2003 pointings with higher counting statistics (observations 4 and 5 in Table 
[1]). Although there is some evidence of intrinsic variability of the thermal component, no 
definite conclusion can be reached at present because of the insufficient statistics. Finally, 
we note that the flux of the MCD component is comparable to that of the PL component 
(see Table ED. 



3.3. Modeling the optical emission 



In order to study the optical emission properties of objects CI and C2 and compare them 
with the HST+VLT photometry, we implemented a model to compute the optical spectrum of 
a binary system with an IMBH taking irradiation effects in t o acco unt. Our calculation relies 
on the same assumptions discussed in ICopperwheat et al.l (120051 ). who recently presented a 
thorough investigation of the infrared-through-optical emission properties of X-ray binaries 
with IMBHs. More speciflcally, we assume that accretion onto the IMBH is fuelled by 
a massive companion fllling its Roche lobe and that the X-ray emission is isotropic; the 
consequences of introducing some degree of beam i ng ar e discussed later on. A standard 
Shakura-Sunyaev disk (e.g. iFrank. King fc Raind 120021 ) is assumed and both the X-ray 
irradiation of the companion (including th e effects of disk shadowing) and the self- irradiation 
of the disk are accounted for. Following ICopperwheat et al.l (120051 ). radiative transfer at 
the donor and disk surfaces is treated assuming a plane-parallel atmosphere in radiative 
equilibrium, i lluminated by t he X-ray flux emitted from the innermost part of the accretion 
disk (see also lWu et al.ll200ll ). In order to keep our treatment simple, the companion star is 
taken to be spherical, neglecting the effects produced by the Roche lobe geometry and also 
those related to the (possible) deformation induced by radiation pressure. Limb and gravity 
darkening were not included. 

The model depends on the masses of the two components, the binary period (which, in 
turn, flxes the orbital separation), the accretion rate and the (unirradiated) temperature of 
the donor, in addition to the inclination angle and the orbital phase. The accretion efficiency 
and the albedo of the surface layers were chosen to be ?7 = GMbh / c^fin = rs/2rin = 0.17 
(where rs is the Schwarzschild radius and T in = 3rs is t he inner disk radius) and a = 0.9 (e.g. 
de Jong, van Paradijs fc Augusteijnj Il996l ) . Following ICopperwheat et al.l (120051 ) , we took 
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the hardness ratio C, = Fx{< 1.5keV)/Fx(> 1.5 keV) = 0.1. The absorption parameters 
in the same two spectral bands were selected as kg = 2.5 and kh = 0.01. The V and B 
magnitudes of the (irradiated) disk plus donor have been computed for several values of the 
parameters of the binary. Each sequence of models, at fixed inclination angle i, accretion 
rate M and donor mass M, corresponds to a track in the color-magnitude diagram (CMD; 
B-V vs. V here) along which the BH mass varies. Only inclination angles smaller than ~ 70° 
are considered because eclipsing effects of the accretion disk on the donor (and viceversa) 
have not been taken into account. On the other hand, from the available data there is no 
positive evidence for eclipses of the X-ray source or the donor in NGC 1313 X-2. The mass 
and luminosity class of the donor fix its (unirradiated) surface temperature Tg//. Different 
tracks have been obtained varying the orbital period Porb which, in turn, determines the 
Roche lobe radius. The maximum allowed period is that for which the Roche lobe radius is 
equal to the donor radius. The computed tracks are compared with the optical (unreddened) 
magnitudes and colors of objects CI and C2 for both the VLT and i/ST observations in § 13.41 

3.4. Objects CI and C2 

In order to constrain the parameters of the binary, we used the optical binary emission 
model introduced in Section O One of the 2003 XMM-Newton pointings of NGC 1313 X-2 
is within 2-3 days from the 1st HST epoch, while another is close to the VLT observation. 
These two XMM-Newton observations are those of 2003 November 25 and December 23; the 
latter was preferred to the observation of December 25 because of the higher statistics. It 
is therefore of interest to compare the V and B magnitudes of objects CI and C2 in these 
two epochs. The variation of the unabsorbed X-ray flux between the same epochs is ~80% 
(see Table H]). At the same time, however, the V and B band magnitudes of CI do not show 
significant evidence of variability. The relative photometric error between the HST and VLT 
data as measured on a sample of field stars is in fact ~ 0.3 mag (see Figure [2]). A similar 
conclusion is reached also for object C2. Again, the magnitude change between the HST 
and VLT epochs is always smaller than the relative photometric error (~ 0.5 mag for objects 
fainter than V=24; see again Figure [2]). 

As an initial guess for the donor parameters in our model we use the values inferred on 
the basis of VLT photometry (M05): an 09-BO V star of ~ 20Mq, T^jf ~ 30000 K for CI, 
and a G-K I star of ~ IOMq, T^ff ~ 4500 K for C2. As mentioned above, this is consistent 
with what inferred from the HST photometry. The donor star in the case of of object CI 
is assumed to be on the zero age main sequence and in contact with the Roche lobe (i.e. 
its radius is equal to the Roche lobe radius). Possible evolutionary effects or disturbances 
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caused by the intense mass transfer are not accounted for. Results for object CI are shown 
in Figure H] for two different values of M, chosen in such a way to match the XMM-Newton 
flux measured in the two observations of 2003 November 25 and December 23. The tracks 
on the CMD diagram are in agreement with the observed V band magnitude and (B-V) 
color of object CI for Porb — 1-7 d, M ~ ISM© and Te// ~ 25000 K (corresponding to an 
early B main sequence star). Taking into account for current uncertainties on both color and 
magnitude, the companion mass and temperature, and the orbital period may vary in the 
ranges 10 < M/M© < 18, 20000 K < T^fj < 30000 K, and 1.5 < Porb < 2 d, respectively. If 
the donor makes contact with its Roche lobe and 0.1 < M/Mbh ^ 0.8, the orbital period 
becomes a function only of the donor radius (or mass). In fact, whe n combining together 
the I II Kepler's law and the expression for the Roche lobe radius (IFrank. King fc Raine 



20021 ) . the dependence on Mbh disappears. As shown in Figure HI the VLT and HST 



photometric points intersect the corresponding tracks at about the same value of the IMBH 
mass. This value increases with increasing inclination angle. Results shown in Figure H] refer 
to orbital phase zero (superior conjunction). The variation in the V (B) band between the 
1st HST epoch and the VLT one is ~ 0.23 (~ 0.25), consistent (within the errors) with what 
observed. Thus, although in these systems X-ray irradiation is very intense, the induced 
optical variability is not very large owing to the high intrinsic emission of the massive B 
donor. The calculation for phases 0.25 and 0.5 gives results similar to those obtained for 
phase 0, typically within 0.15 mag. Thus, this is the expected amplitude of the modulation 
possibly induced by the orbital motion. It is interesting to note that this result is consistent 
with the degree of variability observed in the V band between the two HST observations 
(~ 0.1 mag; see Figure [2]). 

A direct comparison of the three cases illustrated in Figure H] shows that relatively large 
values of the inclination angle {i > 50° — 60°) are required in order to obtain the correct 
optical flux for a BH mass Mbh ~ 12OM0. At lower inclination angles (top and middle 
panels of Figure H]), the black hole masses needed to reproduce the optical magnitudes and 
color of object CI are too small for the X-ray flux to be below the Eddington limit (if the 
emission is isotropic). Therefore, unless the Eddington limit can be circumvented, the binary 
system is expected to have a significant inclination angle. We note that, for any value of 
i, the BH mass inferred from Figure S] is always > 70Mq. This absolute lower limits for 
Mbh is still compatible with a very small beaming (> 0.6). In this case the general picture 
discussed above reasonably continues to hold. On the other hand, if the ULX emission is 
more seriously beamed, the BH mass could be smaller: for a beaming factor ~ 1/6, Mbh can 
be as small as 2OM0 without exceeding the Eddington limit. In this case, we expect no X-ray 
irradiation of both the disk and the companion, being the emission coUimated away from 
the orbital plane. To test this possibility we computed a new sequence of models, following 
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the same procedure outlined above, but switching off the disk/donor irradiation. It turns 
out that it is possible to reproduce the correct magnitude and color, although the donor is 
now less massive and cooler. However, this implies that the star is too small to fill its Roche 
lobe and thus accretion can not proceed through Roche lobe overfiow. Wind accretion may 
still be possible, although it seems unlikely that it can produce the required value of M . 

The situation for object C2 is reversed. We explored the parameter space by varying 
the donor mass and orbital period, but did not find any combination of values which could 
reproduce the data in the framework of isotropic emission. In particular. X-ray irradiation 
causes the (B-V) color to always exceed the observed one. On the other hand, a massive 
and very extended K-type supergiant (M ~ I6M0, Te// ~ 4000 K, Porb ~ 800 days) would 
have properties consistent with those of object C2 if the black hole mass is ~ 2QMq. This 
of course requires a (moderate) beaming. We checked that a beaming factor of ~ 1/6 is 
enough and that the companion fills its Roche lobe. The optical magnitudes are correctly 
reproduced because the (unirradiated) disk contribution becomes negligible in comparison 
with the star intrinsic luminosity. However, in this case practically no variation in the optical 
is expected in response to an increase of the accretion rate, and the predicted magnitudes of 
C2 are constant. This is in contrast with the evidence of variability observed in the V band 
between the two ifST observations (~ 0.1 mag; see Section [2^ and Figure [2]), although some 
variations may be induced also by the orbital ellipsoidal modulation of the donor (which we 
did not take into account). 



4. Discussion 

Although present data do not allow to reach a definite conclusion on the actual coun- 
terpart of the ULX NCG 1313 X-2, some firm points may be derived from the analysis 
presented in the preceeding sections. If CI is the counterpart, as it seems more likely, our 
model indicates that this is an IMBH X-ray binary with a relatively massive main sequence 
donor which fills its Roche lobe. Taking a black hole mass of ~ I2OM0, as required to 
account for the observed X-ray fiux in terms of isotropic emission, the donor mass is in the 
interval 10 — I8M0 (taking photometric uncertainties into account, § 13. 4p . This is larger 
than the maximum main sequence mass of the parent stellar association, ~ 8 — 9M p^ , esti- 



mat ed using multi c olor p hotometry and isochrone fitting by iPakuU. Grise fc MotchI (120061 ) 
and iRamsey et al.l ( 20061 ). However, considering that the lower bound for the donor mass is 
IOMq, the difference is small. We note also that, if CI is the counterpart and C2 belongs 
to the same stellar association, the estimated masses of the two stars correctly places them 
on (or close to) and out of the main sequence, respectively. If the counterpart is C2 then 
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the source is a binary formed by a late type, massive supergiant and a stellar mass black 
hole with beamed X-ray emission. However, this scenario has some shortcomings. First, 
it predicts little if no optical variability, and this is in apparent contrast with the varia- 
tions seen in the two HST observations. Second, the duration of the supergiant phase for a 
^ I5M0 star is very short (a fewx 10^ yr), making the possiblity of catching the binary is 



such an evolutionary stage not very likely (jPatruno fc Zampierill2006l ). Third, the accurate 



astrometry of the field made possible by the cross-identification on the Chandra and HST 
images of a background galaxy appears to strengthen the association of NGC 1313 X-2 with 
object CI, ruling out object C2 (J.-F. Liu, private communication). 

For object CI, an orbital modulation of amphtude ~ 0.15 is expected because 
of orbital inclination and X-ray irradiation effects. This modulation is superimposed to a 
comparable variation caused by changes in the irradiating X-ray flux (~ 0.2 mag). In this 
respect, it is interesting to note that similar variations in the observed B band VLT-I- Subar u 



photometry of object CI have been recently reported also by lPakuU. Grise fc MotchI (120061 ). 
consistent with our findings. In principle, with a sufficient and suitably spaced number 
of observations, the orbital modulation can be singled out and measured with large area 
ground telescopes or HST. The detection of this modulation would lead to the unambiguous 
determination of the orbital period of the binary. This, in turn, would allow us to definitely 
discriminate between CI and C2 and, most importantly, to constrain the mass ratio of NGC 
1313 X-2 and, eventually, the mass of the black hole. 
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Table 1. Observation log of the XMM-Newton EPIC pn pointings and of the 
VLT+FORSl and HST+ACS photometric observations of NGC 1313 X-2. 





Instrument 




Obs. Id. 


Date 


Exposure 


OTP 


Filter 


1 


XMM-Newton EPIC 


pn 


0106860101 


2000-10-17 


31637s 


20600s 


Medium 


2 


XMM-Newton EPIC 


pn 


0150280101 


200.3-11-25 


8365s 


1087s 


Thin 


3 


XMM-Newton EPIC 


pn 


0150280201 


2003-12-09 


5620s 


Os 


Thin 


4 


XMM-Newton EPIC 


pn 


0150280301 


2003-12-21 


10334s 


8272s 


Thin 


5 


XMM-Newton EPIC 


pn 


0150280401 


2003-12-23 


14094s 


3600s 


Thin 


6 


XMM-Newton EPIC 


pn 


0150280501 


2003-12-25 


15282s 


1668s 


Thin 


7 


XMM-Newton EPIC 


pn 


0150280701 


2003-12-27 


16666s 


Os 


Thin 


8 


XMM-Newton EPIC 


pn 


0150280601 


2004-01-08 


14756s 


7696s 


Thin 


9 


XMM-Newton EPIC 


pn 


0150281101 


2004-01-16 


7034s 


4208s 


Thin 




VLT+FORSl 






2003-12-24 


840x2 




B 




VLT+FORSl 






2003-12-25 


600x2 




V 




H5T+ACS(I epoch) 






2003-11-22 


580x2 




F555w 




HST-tACS 






2003-11-22 


630x4 




F435w 




i?5T+ ACS (II epoch) 






2004-02-22 


600x4 




F555w 



''Good Time Intervals in which the total off-source count rate above 10 keV is < 1.0 counts 
s-i. 



Table 2. Spectral analysis of XMM-Newton EPIC pn data of NGC 1313 





Date 


Count rate 


Flux'' 






kTMCD [kcV] 


r 


xl^d/d.o.f. 


1 

2 
4 
5 
6 
8 
9 


2000-10-17 
2003-11-25 
2003-12-21 
2003-12-23 
2003-12-25 
2004^01-08 
2004^01-16 


0.24 
0.67 
0.80 
0.89 
0.56 
0.39 
0.34 


4 00+°-^^ 

r 07+0.28 

9.19i?i° 

4.38ii:;i 


0.93 

0.63 
0.94 
0.87 
0.47 
0.45 


6.52il:g^ 
8.30l°2°^ 

14 qq+5 06 
16.3llJ^2f 

'■'-^-l.SZ 


16+°°2 

15+0-03 

-0.03 
1 5+0-0'' 

n 1 O + 0.03 
'^■'-^-0.02 

17+0-02 


9 o+O.l 
^■•^-0.1 
9 O + 0.2 
^••-'-0.2 

1 9+0-1 
1-^-0.1 

1 S+0-1 
1-8-0.1 
9 1 +0.2 
^-1-0.2 

2 5+0-1 

9 o + O.l 
^•■^-0.1 


1.24/75 
1.06/44 
1.11/144 
0.99/140 
0.99/42 
0.82/67 
1.01/66 



"^Njj has been frozen at 4.02 X 102lcm~2 

I'Unabsorbed flux in units of 10^12 erg em^2 g-i 

"^Ratio of the (unabsorbed) fluxes of the MCD and PL components. 

'lUnabsorbed luminosity, in units of 10^^ erg s~i, for a distance of 3.7 Mpc 
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Table 3. Observed magnitudes and colors of the two candidate optical counterparts of 

NGC 1313 X-2. 





Filter 


CI 


02 


HST I epoch 


B 


23.72±0.04 


26.02±0.04 


VLT 


B 


23.50±0.15 


> 25.2 


HST I epoch 


V 


23.75±0.04 


24.46±0.04 


VLT 


V 


23.60±0.15 


24.10±0.15 


HST II epoch 


V 


23.61±0.04 


24.57±0.04 


HST I epoch 


B-V 


-0.03±0.06 


1.56±0.06 


VLT 


B-V 


-0.1±0.2 


> 11 
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Fig. 1.— HST+ ACS images of NGC 1313 X-2: left F435W (B) band, right F555W (V) 
band. The Chandra error box, the candidate optical counterparts CI and C2 and the field 
sources A and D are shown (following Z04 and M05). 
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Fig. 2. — Top and middle panels: AV = V — V{HST I epoch) and AB = B — 
B{HST I epoch) for objects CI and C2 for the available epochs. The error bars corre- 
spond to 0.3 and 0.5 mag for V and B respectively (see text for details). Lower panel: 
XMM-Newton count rates of NGC 1313 X-2 in the [0.2-10.0] keV range. 
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Fig. 3. — The X-ray lightcurve of NGC 1313 X-2. Fluxes are unabsorbed and refer to the 
[0.2-10.0] keV energy interval (see Table [2]). Einstein, ROSAT and ASCA points are taken 
from Z04. The insert refers to the more recent XMM-Newton data (Table [2]) and the arrows 
mark the time of the HST and VLT observations (Tabled]). 
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Fig. 4. — Color-magnitude diagram for the (irradiated) disk plus donor model for Port — 
1.7 d, M ~ ISM© and Te// ~ 25000 K (object CI). Each panel refers to a different 
inclination angle i. The two tracks correspond to M = 3.5 and 6.7 MEdd- These values are 
chosen in such a way to match the XMM-Newton flux measured in the two observations 
of 2003 November 25 and December 23, that are quasi- simultaneous to the HST I epoch 
and VLT observations, respectively. The labels indicate the BH mass. The V magnitude 
and B-V color as obtained from the VLT and the HST observations are also shown [open 
squares). 



